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We report the generation of large coherent vortices via inverse energy cascade in Faraday wave 
driven turbulence. The motion of floaters in the Faraday waves is three dimensional but its horizon- 
tal velocity fluctuations show unexpected similarity with two-dimensional turbulence. The inverse 
cascade is detected by measuring frequency spectra of the Lagrangian velocity, and it is confirmed 
by computing the third moment of the horizontal velocity fluctuations. This is observed in deep 
water in a broad range of wavelengths. The results broaden the scope of recent findings on Faraday 
waves in thin layers [A. von Kameke, et al, Phys. Rev. Lett. 107, 074502 (2011)]. 

PACS numbers: 47.35.Pq, 47.27.Cn, 47.52.+j 



Ripples appearing on the surface of a vertically vi- 
brated liquid are known as Faraday waves. They are 
widely used to study a broad range of phenomena, such 
as pattern formation [IHi], solitons [SUTOl. walkers, or 
bouncing orbiting droplets [11] , extreme wave events [12] , 
and others. It has recently been found [TS that the hori- 
zontal motion of floating particles on the surface of Fara- 
day waves for thin layers of liquids shows several proper- 
ties consistent with the fluid motion in two-dimensional 
(2D) turbulence. This turbulent horizontal transport was 
referred to as Faraday flow. In particular, a Kolmogorov 
^-5/3 spectrum has been reported as well as the Richard- 
son law of pair dispersion [14 . Even though the layer was 
argued to be thin (i.e. smaller than the Faraday wave- 
length), the Faraday flows explored in these experiments 
are indisputably three dimensional (3D) and locally di- 
vergent flows. Some of the 2D turbulence theory assump- 
tions are thus clearly violated and one should have few 
expectations for the emergence of the inverse energy cas- 
cade in this context. 

However it has recently been observed in electromag- 
netically driven turbulence that the behavior of a 3D flow 
can be rendered similar to 2D turbulence in the presence 
of a large scale vortex [15 . Moreover, recent numer- 
ical predictions emphasize that an inverse energy cas- 
cade might be a general principle existing in all turbulent 
flows in both 2D and 3D ^6j. The emergence of such a 
phenomenology in Faraday waves would broaden the ap- 
plicability of features common to 2D turbulent flows to 
surface wave phenomena. But the key ingredients sup- 
porting this proposition are to date still unclear. Indeed 
many features of Faraday flows are still unknown, such 
as the spectral localness of the energy injection in the 
surface ripples, the limitation to thin layers of fluid, or 
the dependence of the floaters motion on the frequency 
of vertical vibration. Besides it would be interesting to 
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discover whether it is possible to generate large coher- 
ent structures, or spectral condensates, in the Faraday 
waves such as those observed in bounded electromagnet- 
ically driven turbulence [TTl [18] . 

In this paper we report new experimental results which 
show that the particle motion on the undulating surface is 
similar to the fluid motion in 2D turbulence. It supports 
the inverse energy cascade or the spectral energy transfer 
from smaller to larger scales. The water depth in these 
experiments is not restricted to thin layers, but rather 
substantially exceeds the Faraday wavelength. Moreover 
the vertical acceleration ranges from the Faraday insta- 
bility threshold up to the droplet nucleation threshold 
where the ripples are a couple of millimeters high. Such 
a configuration rules out any 2D assumption on the fluid 
motion. If the flow is bounded, the inverse cascade leads 
to the accumulation of energy at the boundary box scale 
and the generation of a large vortex coherent over the 
whole domain. 

Experiments are carried out in a circular container (178 
mm diameter, 30 mm deep) as shown in Figure [l] The 
container is filled with a liquid whose depth is larger than 
the wavelength of the perturbations at the surface. The 
monochromatic forcing is varied in the range of frequen- 
cies /o = 30 — 60 Hz which results in Faraday wavelengths 
A in the range A ~ 8 — 14 mm. The dominant frequency 
of the excited surface ripples is at the first subharmonic 
of the excitation frequency, / = /o/2. We define the 
supercriticality as a = (a — athj/dth , where a is the am- 
plitude of the vertical acceleration imposed by the shaker 
and ath is the threshold for the parametric generation of 
Faraday waves. 

A diffusing light imaging technique is used to visual- 
ize simultaneously the surface ripples and floating tracer 
particles, as seen in Fig. |2|a). A 2% milk solution in 
water provides sufficient contrast for the parametrically 
excited waves to be observed and allows high resolution 
movies (800x800 pixel) to be recorded in the range of 
speeds 80 - 120 fps using the Andor Neo sCMOS camera. 
Black floating particles are added on the fluid surface to 



2 



iflv camera 



side view f 

1 1 1 1 1 1 1 1 1 



top view 




fluid container with 
transparent bottom 



LED light 
source 



electro-magnetic 
shaker 



FIG. 1. Experimental setup for the surface wave driven 
turbulence. 

visualize the horizontal motion of the liquid. Particles are 
made of carbon glass and have been plasma treated to 
reduce their intrinsic hydrophobicity. The use of surfac- 
tant and plasma treatment ensures that particles do not 
aggregate at the surface. This is illustrated in Fig. |2|a) 
with particles homogeneously distributed at the surface 
of the wave field. 

It has recently been shown that the Faraday waves 
properties are consistent with those of an ensemble of 
bounded oscillating solitons, or oscillons [10]. Figure[2|b) 
shows examples of the oscillons' trajectories to be com- 
pared with a few particle tracks measured using the same 
tracking code. The horizontal motion of oscillons is ran- 
dom in time but confined to a disordered lattice as seen 
from the black trajectories shown in Fig. [2|b) [20 . The 
particles also wander erratically, but their trajectories 
show substantially larger excursions in comparison to os- 
cillons. A particle usually visits several oscillon sites 
which randomize its trajectory. Although the particle 
and oscillons trajectories are characterized by completely 
different length scale, the oscillon mobility is an essential 
feature of Faraday flows. In the case of a perfect oscil- 
lonic crystal, in which the relative position of the nodes is 
frozen, the horizontal mobility of particles is zero, show- 
ing that the horizontal particle transport is strongly cor- 
related to the oscillon motion. Such crystals are obtained 
by adding minute amounts of bovine serum to water [10] . 

Though the oscillon lattice in Fig. [2|a) seems disor- 
dered, the wave number spectrum of the ripples is very 
narrow, as seen in Fig. |2|c). This is consistent with the 
fact that oscillon motion is restricted in space to about 
half of the lattice characteristic wavelength [20 . To gen- 
erate 2D turbulence one needs to inject energy into the 
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FIG. 2. (a) Diffusive light image of the fluid surface elevation 
including floating particles. The particles diameter here is 300 
microns 19 , the peaks and troughs of the oscillons appear as 
dark and white blobs, (b) Trajectories of oscillons (black) 
having a phase maximum, and trajectories of 8 particles (col- 
ored) tracked for 2 s (200 periods of the shaker), (c) Wave 
number spectra of the surface elevation at different vertical 
accelerations (a = 0.7 g, 1.2 g, and 1.6 g) and /o = 60Hz. (d) 
Frequency power spectra of the Lagrangian velocity (averaged 
over ^ 1000 trajectories) for different vertical accelerations at 
/o = 60Hz. 

horizontal fluid motion at some intermediate range of 
scales in a localized wave number domain. The horizontal 
mobility of oscillons is the manifestation of a momentum 
transfer along the water surface; this transfer is restricted 
to a narrow A:-domain. 

Fig. |2|d) shows frequency power spectra of velocity of 
particles moving along their Lagrangian trajectories. The 
spectra show peaks close to the Faraday frequency (30 
Hz, the first subharmonic of the excitation frequency), 
and a broad low-frequency band whose energy grows with 
the increase in vertical acceleration. Since energy is in- 
jected at / = 30 Hz, the growth and the broadening 
of the low-frequency band indicates a nonlinear energy 
transfer to slower temporal (larger spatial) scales. This 
transfer, as will be shown below, is due to the inverse 
energy cascade. 

Fig. [sj^a) shows many trajectories of particles on a 
densely seeded fluid surface, the presence of multiple 
length scales in this Faraday flow is visible. In the follow- 
ing. Particle Image Velocimetry (PIV) is used to charac- 
terize the Eulerian velocity field [21 . 

Figure |3jb) shows wave number spectra of the horizon- 
tal kinetic energy measured at /o = 60 Hz at three ver- 
tical acceleration levels, a = 0.7g, 1.2g, and 1.6g, corre- 



3 




FIG. 3. (a) Particle trajectories followed over 5 periods of the shaker oscillations at a = 1.6g and /o = 60Hz in a 8 x 8 cm^ 
field of view, (b) Wave number spectra of the horizontal kinetic energy at different levels of the vertical acceleration (a = 0.7 
g, 1.2 g, and 1.6 g) and /o = 60Hz. (c) Kinetic energy of the flow Eo versus vertical acceleration a. Grey bar indicates the 
threshold of parametric instability. 



spending to supercriticalities in the range a = [0.17, 1.7] 
{ath = 0.6g at 60 Hz). The spectrum is close to the 
Kolmogorov-Kraichnan theory power law of Ey(k) ex 
k~^^^ at wave numbers k < 1500 m~^. At higher wave 
numbers, k > 1500 m~^, spectra are somewhat steeper 
than expected from the direct enstrophy cascade fit of 
probably due to higher damping at large wave num- 
bers. If we define the turbulence forcing wave number 
from the position of the kink on the turbulence spectrum, 
e.g. kf ^1550 at 60 Hz, Fig. [sj^b), it appears that 
the turbulence forcing wave number is roughly twice that 
of the surface elevation wave number (Fig.j^Jc)). The en- 
ergy injection scale for the horizontal particle transport 
is thus related to the oscillon size which is about half of 
a period of the Faraday wave. 

The kinetic energy of the horizontal flow can be es- 
timated as Eq = Ey{k)dk. Here, kiow is the low- 
est wave number determined by the field of view. In 
Fig. (sjc), £^0 is shown as a function of the vertical ac- 
celeration a. By changing the vertical acceleration from 
0.7g to 2.4g, we can vary the kinetic energy in the flow 
by over two orders of magnitude. 

Similar measurements were performed at several fre- 
quencies of the vertical vibrations: 30 Hz, 45 Hz, 60 Hz, 
Fig. [i^a). A k~^^^ power law is consistently observed in 
the energy spectra at lower wave numbers. A distinct 
kink is present in all spectra, and it is associated with 
the forcing wave number. The turbulence forcing wave 
number kf decreases with the decrease in frequency, in 
accordance with the capillary-gravity wave dispersion re- 
lation. A k~^^^ range in the horizontal velocity spectra 
is consistent with the Kolmogorov-Kraichnan theory. 

An important feature of 2D turbulence is the pres- 
ence of the inverse energy cascade, or the transfer of en- 
ergy from smaller to larger scales in the inertial range 
Vf < r < riow A direct way to prove the existence of 
the inverse energy cascade is to compute the third mo- 
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FIG. 4. (a) Energy spectrum in the surface wave driven 
turbulence at different frequencies (30 Hz, 45 Hz, 60 Hz), and 
at similar supercriticalities. Colored vertical bars show the 
forcing ranges, (b) Third moment of the longitudinal velocity 
increments measured at /o = 30 Hz at different supercritical- 
ities a = l(a =0.4g), and a — 2(a =0.6g). 

ment of the velocity increments across a distance r in 
a flow, Sv{r) = v{l -\- r) — v{l). The third order mo- 
ment SsL =< > (where index L refers to the 
longitudinal velocity component with respect to r, and 
angular brackets denote averaging over all possible po- 
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FIG. 5. (a) Particles trajectories followed over 4 periods 
of the shaker oscillations in a spectrally condensed flow in a 
square box. The large scale vortex is 10 cm in diameter. The 
particles diameter is 50 microns, (b) Kinetic energy spectrum 
in the spectrally condensed turbulence. 

sitions / within the flow field, and in time, over many 
realizations) is related to the spectral energy flux e via 
the Kolmogorov flux relation. In the inertial range of 2D 
turbulence this relation reads as [T71 [22] 

SzL = lev. (1) 

Fig. Qb) shows third order structure functions com- 
puted from the velocity fields measured in the surface 
wave driven turbulence at /o = 30i72: for two different 
vertical accelerations. In both cases S^l is a positive 
linear function of r. The spectral energy flux e can 
be computed from Eq.(l). The derived values are in 
agreement with the spectral power densities derived 
from the spectra of Fig. [4]^a). 

The existence of the inverse energy cascade in a 
bounded flow is a prerequisite for the accumulation of 
the turbulence energy at the scale of the boundary size 
pS'. But the spectral condensation, or the generation of 
large coherent structures, is dependant on the boundary 
size, the energy flux, or the damping rate. It has been 
demonstrated in experiments in electromagnetically 



driven turbulence (EMT) that in a square box, whose 
size is smaller than the dissipation scale, turbulent 
inverse energy cascade can drive large scale vortices 
which are coherent over the entire domain To 
test if the inverse energy cascade can generate coherent 
vortices in the surface wave driven turbulence, we 
produce Faraday flows in a square container (110 x 110 
mm^) at the excitation frequency of /o ~ 30 Hz, which 
gives the forcing scale of If ~ 7.7 mm. Fig. |5|a) shows 
trajectories of particles in this smaller square container. 
A large vortex occupying the entire box develops in the 
square container about 60 seconds after switching on 
the shaker. The kinetic energy spectrum of spectrally 
condensed flow is substantially steeper than the Kol- 
mogorov spectra in Fig. [sj^b), E{k) ex as shown in 
Fig. [sjb), similarly to the observations in EMT [iT] , 

To conclude, experimental results support the exis- 
tence of the inverse energy cascade fueled by Faraday 
waves. Faraday flows represent thus a new versatile tool 
of laboratory modeling of 2D turbulence. It is remark- 
able that the inverse energy cascade, initially thought of 
as a process intrinsic only to idealized 2D turbulence is 
found in a variety of physical systems where its existence 
could not possibly be expected a priori [I3j [T5j [24] . In 
3D turbulence, the inverse energy cascade was recently 
found to naturally coexist with the direct energy cascade 
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